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ABSTRACT 


The response of a thin, rigid viscoplastic plate subjected to a spatially 
axisymmetric Gaussian ideal impulse loading was studied analytically. The 
Gaussian ideal impulse distribution instantaneously imparts a Gaussian initial 
velocity distribution to the plate, except at the fixed boundary. The plate 
deforms with monotonically increasing deflections until the initial dynamic 
energy is completely dissipated in plastic work. The simply supported plate of 
uniform thickness obeys the von Mises yield criterion and a generalized consti- 
tutive equation for rigid, viscoplastic materials. For the small deflection 
bending response of the plate, neglecting the transverse shear stress in the 
yield condition and rotary inertia in the equations of dynamic equilibrium, the 
governing system of equations is essentially nonlinear, A proportional loading 
technique, known to give excellent approximations of the exact solution for the 
uniform load case, was used to linearize the problem and obtain analytical 
solutions, in the form of eigenvalue expansions. The linearized governing equa- 
tion required the knowledge of the collapse load of the corresponding static 
problem. 

The effects of load concentration and an order of magnitude change in the 
viscosity of the plate material were examined while holding the total impulse 
constant. In general, as the load became more concentrated, the peak central 
velocity increased and the time for plate motion to cease increased. For the 
less viscous plate, these increases of velocity and time were more pronounced. 
The final plate profile became more conical as the load concentration increased, 
but did not approach the purely conical shape predicted for the point impulse by 
the rigid, perfectly plastic analysis with the Tresca yield criteria. Profiles 
of the less viscous plate were influenced more by the load concentration. 
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deviator stress tensor at initial yield 
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time for motion to cease 
dynamic component of velocity 
steady component of velocity 
initial velocity 
nondimens ional plate velocity 
transverse deflection of the plate 
transverse coordinate 


plate geometry and material constant, sec 


nondimens ional Gaussian shape parameter 
material constants 

harmonic and biharmonic operators in cylindrical coordinates 
final center deflection 
strain rate tensor 

eigenvalues determined from equation (A7) 
mass density per unit area of the plate 
nondimens ional radial coordinate 
stress tensor 

yield stress in simple tension 
function defined by equation (3) 
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function defined by equation (A12) 
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INTRODUCTION 


This paper presents the results of an analysis of the small-deflection 
bending response of a simply supported circular plate of rigid, viscoplastic 
material subjected to a spatially axisymmetric Gaussian ideal impulse. The 
effects of load concentration and an order of magnitude change in the viscosit 
of the plate material are examined while holding the total impulse constant. 
Approximate expressions are developed for the time at which plate motion 
ceases, the final shape of the plate, and the final central displacement. 

Although there have been a number of papers (refs. 1, 2, 3) which permit 
a time variation of the load, there have been few papers which consider a 
radial variation other than linear (refs. 3, 4). The only general spatial 
distribution of load which has received significant analytical attention is tt 
Gaussian distribution. By varying a single parameter, this general distribu- 
tion can span the extremes from the point load to the uniformly distributed 
load. This versatility was recognized by Sneddon (ref. 5) who approximated 
the dynamic loading of a projectile on a thin, infinite elastic plate by a 
Gaussian distribution of pressure. Madden (ref. 6), in his study of shieldinj 
of space vehicle structures against meteoroid penetration, related the 
meteoroid-shield debris loading of the main vehicle wall to a Gaussian initial 
velocity distribution. The first study of this loading on a plastic plate was 
by Thomson (ref. 7). He obtained the solution of a rigid, perfectly plastic 
plate of material obeying the Tresca yield condition subjected to an initial 
impulse of Gaussian distribution. Weidman (ref. 2), in/considering the 
response of simply supported circular plastic plates to distributed time- 
varying loadings, presented an example case of a radial Gaussian distribution 
of pressure with an exponential decay. The plate material was also rigid, 
perfectly plastic obeying the Tresca yield conditions. 

A generalized constitutive equation for rigid, viscoplastic materials is 
presented in the next section. Material elasticity is neglected in order to 
simplify the analysis, as is frequently done in theoretical investigations of 
dynamic plastic response of structures. Rigid-plastic analyses are generally 
believed to be valid when the dynamic energy is considerably larger than the 
maximum energy which could be absorbed in a wholly elastic manner and the 
duration of loading is short compared with the fundamental period of vibratio 


LINEARIZATION OF THE GENERALIZED CONSTITUTIVE EQUATIONS 
FOR RIGID, VISCOPLASTIC MATERIALS 


Perzyna (ref. 8) developed a generalized constitutive equation for rate 
sensitive plastic materials by incorporating a general function in the 
relationship to take the place of the yield function as used by previous 
researchers (Hohenemser and Prager, ref. 9; and Prager, ref. 10). Utilizing 

the definition of the second invariant of the stress deviator, ^ 

the yield function is expressed as 
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( 1 ) 


F = 



where S-jj is the stress deviator tensor and k is the yield stress. The 
generalized constitutive equation proposed by Perzyna is 


3F 


£ ii - 


ij 


where £.. is the strain rate tensor, 

$(F) =0 if F < 0 
#(F) t 0 if F > 0 

and y° denotes a physical constant of the material. 


( 2 ) 


(3) 


Perzyna (ref. 11) has shown that the generalized constitutive equation 
for viscoplastic materials reduces to the constitutive equations of an incom- 
pressible, perfectly plastic material first considered by von Mises and to the 
flow law of perfect plasticity theory. As in the theory of perfectly plastic 
solids, convexity of the subsequent dynamic loading surfaces and orthogonality 
of th e inelastic strain-rate vector to the yield surface follow from Drucker's 
postulates defining a stable, inelastic material with inclusion of time- 
dependent terms (Perzyna, ref. 8). 


A method of linearizing boundary-value problems in the theory of visco- 
plastic solids is described by Wierzbicki in reference 12. In this method, as 
shown graphically in figure 1, the concept of proportional loading is used to 
relate the state of stress S-ji on the initial yield surface F = 0 to sub- 
sequent states of stress, namely, proportional loading requires the direction 
cosine tensor of the state of stress in deviatoric space to be independent of 
time: 



(4) 


This is a reasonable approximation for axisymmetrically loaded simply supported 
circular plates because the plate -center and boundary are automatically pro- 
portionally loaded, that is, the bending moments must always be equal at the 
plate center and the circumferential bending moment must always be zero at the 
plate boundary. 

Utilizing equation (4) , the generalized constitutive equation (2) 

>ecomes 
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where the viscosity constant y = Y°/2k. For this analysis,, the linear form 

$(F) = F (6 


is chosen. This simplified constitutive equation still is nonlinear in stress' 
However, in the solution of dynamical plate and rotationally symmetric shell 
problems, the constitutive equation (5) with the linear function $(F) - F 
produces full linearization of the governing equations. 

For the problem of a uniformly loaded, simply supported circular plate 
with $(F) = F, Wierzbicki (ref. 12) has shown that the approximate solution 
obtained using the proportional loading hypothesis agrees very well with a 
numerical finite-difference solution of the exact equations. The solution of 
the linearized problem also agrees well with experimental data on impulsively 
loaded plates by Florence (ref. 13). 

For the linear function equation (5) becomes 


hi ■ i (s u - V 


(7 


where equation (7) is really a flow relation for a given structure rather 
than a constitutive equation describing a given material (ref. 14). 


GOVERNING EQUATIONS, BOUNDARY AND INITIAL CONDITIONS 


A Gaussian ideal impulse is suddenly applied to the entire surface of a 
rigid, viscoplastic plate of radius R and thickness 2h resulting in an 
initial velocity distribution described by 


T 2 2 

TI/ I -a r 

V (r , 0) = ~ e 


(8 


where I is the impulse per unit area at the center of the plate and y is 
the mass density per unit area of the plate middle surface. The boundary of 
the plate at r = R is simply supported. The geometry of the plate and 
initial velocity are shown in figure 2. 

The parameter a in the distribution function is a shape parameter whicf 
controls the distribution of the impulse. For a = 0 equation (8) describes 
a uniform impulse; and as a I -* 00 equation (8) describes a point 

impulse at the plate center. 
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The internal forces and moments acting on a typical plate element are 
shown in figure 3. If rotary inertia is neglected, but transverse inertia 
taken into account, the equations of motion are 


fe ™ - - $ 

(9) 

h (rM r>. - M * - r< ! 


Utilizing the Lov e-Kir cho f f hypotheses, the curvature-rate-moment rela- 
tions, derived from the linearized constitutive equation, equation (7), are 


K - t(* r - V ' <*. - V 1 
** = jr [(a** - « r ) - ( 28 * - M r )l 


( 10 ) 


where B = /3 y/2h. M r and are moments satisfying for any r the 

equation of the initial yield surface 


M 2 - M Mr + M 2 = M 2 (11) 

r r 4> <p o v 

M Q = 0 Q h^ is the yield moment of the plate material and 0 Q is the yield 
stress in simple tension. 

• . 

For small deflections of the plate the curvature rates K r and are 

related to the deflection rate w by 
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( 12 ) 


Equations (9), (10), and (12) form a linear parabolic system of partial 
differential equations with six unknown functions — Mj., M(j>, Q, w, 1^, and 
— plus the unknown static moment distribution M r and M^. 

By eliminating all unknowns except w, the system of governing equations 
can be reduced to the single, fourth-order equation 


where 



(13) 
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The right-hand side of equation (13) represents the internal force distribu- 
tion at the initiation of collapse in the static case. 

Let pi denote the static load-carrying capacity of the plate, then the 

„ a 2 r 2 

right-hand side of equation (13) can be replaced by -p' e and the 

governing equation becomes 
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This method of solution, proposed by Wierzbicki (ref. 12), has the impor- 
tant property of replacing the unknown static moment distribution M r and M<j>, 
whose explicit formulas are not known for the von Mises yield condition, by 
the static load-carrying capacity p^. Thus, the need for explicit formulas 
has been reduced to finding the value of a constant, Pq, corresponding to a 
particular value of the shape parameter, a. The determination of the load- 
carrying capacity, p^, of a circular plate under a Gaussian distribution of 
pressure is presented in reference 15. 

Define the dimensionless quantities 
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2 3 ^ 1 3 

and let V - — - + — Then the final form of the governing equation, 

3p 2 P 3P 
equation (14) , is 


4 . 3 9v 0 /o -c, -3p“ 

V v+^ag— = -2/3F' e 


(16) 


where a = yBR /M q . 

The boundary conditions of the simply supported plate are 

m = n, q = 0 at p = 0 


(17) 


v = m = 0 at p = 1 
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Using equations (10), (12), and (9), equations (17), in terms of 
rate of deflection become 


lim 
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(18) 


For the Gaussian ideal impulsive loading the plate is initially flat and the 
initial velocity has a Gaussian distribution 


w(p, 0) = 0; v(p, 0) 


II e v 
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(19) 


where I’ 



RESULTS AND DISCUSSION 


The solution to the governing equation, equation (16) , with associated 
boundary and initial conditions, equations (18) and (19) are presented in 
the Appendix. The effects of load distribution and plate viscosity on plate 
response are examined in this section while holding the total impulse constant. 

2 

The impulse amplitude, I' = IR /M 0 , sec, at the plate center is related 
to the total impulse, iT, and distribution parameter, $, by the relation 


I* 



( 20 ) 


T 

I 1 -3 

For comparison purposes the total impulse is held constant at = 1 x 10 

o 

sec. The impulse becomes more concentrated at the center of the plate as 8 
Ls increased and the amplitude grows almost linearly as 8 becomes large. 

For 8=0, the impulse has a uniform distribution. 

The graphical results were obtained by programing the solution (equations 
(A14) and (A15)) and summing the series term-by-term. The rapidly convergent 
series with l/A^ and l/A^ factors did not present any computational diffi- 
sulties; however, the last series in the velocity expression equation (A14) 
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has a lA n factor and prohibited the calculation of velocity-time histories 
for small 3 and t. For t =0 the series is slowly convergent. 

A representative plot of the plate central velocity is shown in figure 4 
for 3 = 10 and viscosity parameter a = 1 X 10” 2 see. The initial central 
velocity is seen to rapidly decline during the first 0.025 msec after which 
the velocity more slowly tends to zero. 

The plate is seen to deform monotonically with increasing deflection until 
the initial dynamic energy is completely dissipated in plastic work and the 
plate comes to rest. The deformed profiles of the plate at rest are shown in 
figure 5 for two values of a (1 x 10 , 1 x IQ” 3 ) and various values of 3- 

The profile becomes more conical as the impulse becomes more concentrated and 
the profiles of the less viscous plate (a = 1 x 10” 2 sec) exhibit a wider 
variation, thus are influenced more by the shape parameter 3 than are those 
for the a = 1 x 10” 3 sec case. 


Approximations 

An approximation to the deflection of the plate is obtained from equation 
(A15) by retaining only the first terms of series and using the approximation 


n / r> 2 s n 


c + c d 2 + i- e' 6p2 + (i + a 2 ) y .(-b ffph 
C 1 + C 2 P + 2B P 2 A (2n) n! 


X 1 


The result is 
1 


w(p,t) 


BR 


/3 F' 

c, + c n p 2 ■ 

4- ~ e + 

(i + 

43 

l T 

23 


ff 

2X 0 i 

r ^ 

a 

L 

3a t \ 

3a , /3 I* 


j' 

i 4 + 4F> i 

A 2X i J 

- t 




J 


,n„ 2 x n‘ 


n=l 


(2n) 


2 N n”l 
P ) 

n! 


( 21 ) 


An approximate expression for the time for motion to cease can be obtained by 
setting the derivative of the approximate displacement expression to zero, 

3wl - 0, is 
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Equation (22) is plotted in rigure 6 for 0 <_ 8 <_ 100 and several values of 
a. Equation (22) is an implicit function of & since I* and F* vary with 
3. The effect of 3 diminishes after an initial rapid rise of tf with 
increasing 3. The symboled points represent computed times using the complete 
equation for the velocity, equation (A14) . Equation (22) is a very good 
approximation for the case a = 1 x 10“3 sec. However, except for small values 
of 3, the approximation is poor for the a = 1 x 10" 2 sec case. 

For a -*■ °° equation (22) limits to 

/3 I* 

4F' 

(23) 
I 


and represents the rigid, perfectly plastic case (y °°) with the von Mises 
yield condition. Equation (23) has the same form as Wang's (ref. 16) result 
for the uniform ideal impulse problem using the Tresca yield condition for a 
rigid perfectly plastic material. However, equation (23) gives slightly 
smaller values of tf since p,!, = 6.51 for the von Mises yield condition 
rather than 6 in the ease of the Tresca yield condition. 

The curve labeled Tresca, r.p.p. was obtained from the results of refer- 
ence 7 where a simply supported circular plate of rigid, perfectly plastic 
material obeying the Tresca yield condition and associated flow rule was 
analyzed for a general Gaussian ideal impulse loading. For small 3 the two 
curves differ only slightly, but as 3 grows larger and the impulse becomes 
more concentrated, the two analyses predict drastically different times for the 
plate motion to cease. The Tresca yield condition predcits very large times 
for plate motion to cease, whereas the von Mises yield condition predicts more 
realistic times for concentrated loads. 

The substitution of equation (22) for tf into equation (21) provides 
an approximate expression for the final plate displacements: 



-“p v(p.t) 

BR 



and for the final center displacement 6(0, tf) = w(0,tf): 
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Equation (25) is plotted as a function of 3 for the two values of a in 
figure 7. The approximations are in excellent agreement with the points 
computed from the exact equations for both a = 1 x 10"^ sec and 1 x 10“ ^ sec, 
even though the tf-approximations for the larger a were poor for large 3 
as shown in figure 6. The nondimensional central displacements are shown 
smaller for the a = 1 x 10“^ sec case when, in reality the real displacements 
are larger than for the a = 1 x 10“ 3 sec case. This is caused by a being 
in the denominator of the expression for the nondimensional central displacement. 


Profiles obtained from the approximation, equation (24) , were compared 
with profiles obtained from the exact equation. For a = 1 x 10“ ^ sec, the 
differences between the approximate and exact profiles were negligibly small 
for the entire range of 3 considered, 10“ ^ to 10,000. However, for the less 
viscous plates, a = 1 x 10“^ sec, the differences were not negligible and the 
approximation, equation (23), should therefore be restricted accordingly. 


CONCLUDING REMARKS 


A thin, simply supported rigid, viscoplastic plate subjected to a Gaussian 
ideal impulse has been analyzed within the realm of small deflection bending 
theory. The plate material obeys the von Mises yield criteria and constitutive 
equations due to Perzyna (ref, 11). These considerations lead, essentially, to 
nonlinear equations governing the dynamic response of the thin plate. A pro- 
portional loading hypothesis, proposed by Wierzbicki (ref. 12) and shown to be 
an excellent approximation of the exact solution for the uniform load case, 
was used to linearize the problem and obtain analytical solutions in the form 
of eigenvalue expansions. The linearized governing equation on the velocity 
of the plate required the knowledge of the collapse load of the corresponding 
static problem, that is, the collapse load for the specific load distribution 
parameter, 3* 

The effects of impulse concentration and an order of magnitude change in 
the viscosity of the plate material were examined while holding the total 
impulse constant. In general, as the impulse became more concentrated, the 
peak central velocity increased and the time for plate motion to cease 
increased. For the less viscous plate material, these increases of velocity 
and time, tf, for plate motion to cease are more pronounced. The final plate 
profile became more conical as the load concentration increased, but did not 
approach the purely conical shape predicted by the rigid, perfectly plastic 
analysis with the Tresca yield condition for a point impulse. As the viscos- 
ity of the plate decreases, the shape parameter has more effect on the final 
deformed plate profiles. 
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Approximate expressions were developed for the time at which plate motion 
ceases, tf, the final shape of the plate, and the final central displacement. 
Comparisons with the series solution indicated that the approximations were 
excellent for the a - 1 x 10"^ sec case. The approximation for the final 
central deflection was good for the entire range of shape parameter $.» the 
other approximations were limited in usefulness. 
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APPENDIX 


SOLUTION OF EQUATION (16) BY EIGENVALUE EXPANSION* 


Since the right-hand side of equation (16) is not a function of time, it 
can be solved by means of an eigenvalue expansion method. Substitution of 

v(p,t) = u(p,t) + U(p) (Al) 


into equation (16) results in 

V 4 u(p,t) + | a + V 4 U(p) = - 2/3 F' e”^ P 

which separates into 

V 4 u + | a || = 0 (A2) 

and 

V 4 U = - 2/3 F' e~ ep (A3) 

Equation (A3) is the same as equation (16) except for the absence of the 
inertia term. Thus, U(p) is an equilibrium solution of equation (16) with the 
same boundary conditions, equations (17). 1 

The solution to equation (A3) satisfying the boundary conditions, equa- 
tions (18), is 
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Eor more details, the reader can consult "Gaussian Impulsive Loading of 
Rigid Viscoplastic Plates," by R. J. Hayduk, Ph. D. Thesis, Virginia Polytech- 
nic Institute and State University, Blacksburg, Virginia, 1972. 
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A general solution due to Wierzbicki (ref. 12) satisfying equation (A2) 
and all prescribed boundary conditions can be written in the form 


T -(2A /3a)t 

u(P,t) . I A [^ttJJ^p) - ra„)I„W„P)] e n 
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n o n r 
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(A6) 


where J Q (x) and I Q (x) denote the Bessel functions of the first kind of real 
and imaginary arguments. The solution (A6) identically satisfies boundary con- 
ditions (18 a, b, and d). The eigenvalues, A n , are roots of the fo llowing 
transcendental equation stemming from the boundary condition (18c) of zero 
bending moment at the plate edge 


W J 1 ( V + WV V - ‘•WWV - 0 


(A7) 


The only remaining unknowns in the solution are the series coefficients 
These coefficients are evaluated from the initial condition (19), that 

is, 

-sp 2 


v(p,0) = u (p , 0) + U(p) = 
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and after substituting equation (A6) fof u(p,0) there results 


T A 1 [I (A )J (A p) - J (A )I (A p)] = * U (p) + ~ e 3p 
n o n o n K o n o n a 


(A9) 


n=l 


The coefficients A^ can be determined from (A9) by virtue of the orthogonal- 
ity of the system [I Q (A n ) J 0 (A n p) - J Q (A n ) I 0 (A n p) ] on the interval [0,1] where 
p is used as a weighting function. Therefore, coefficients A^ can be 
determined as 

1 T 1 -Rn 2 

n (pU(p) - p — e PP )[I (A )J (A p) - J (A )I (A p)] dp 
0 _ _ a o n o n o__ n o n (A10 

J. p[I (A )J (A p) - J (A )I (A p)] 2 dp 

J 0 o n o n o n o n 

where U(p) is defined by equation (A4) . The resulting coefficients are 
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When equations (A4) and (A6) are summed and equation (All) is used, the 
complete solution becomes 


v(p,t) 
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The displacement of the plate is determined by integrating (A14) with 
respect to time. Taking the initial condition of zero displacement into 
account, the displacement becomes 
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Equations (A14) and (A15) represent the complete solution for the velocity 
and displacement of the plate. In the limit as 8 0, the Gaussian ideal 

impulse becomes the uniform ideal impulse and this solution reduces to the 
solution presented by Wierzbicki (ref. 12). 


REFERENCES 


1. Perzyna, P.: Dynamic Load Carrying Capacity of a Circular Plate. Arch. 

Mech. Stos., Vol. 10, No. 5, 1958, pp„ 635-647. 

2. Weidman, D. J.: Response of a Plastic Circular Plate to a Distributed 

Time-Varying Loading. Ph. D, Thesis, Virginia Polytechnic Institute 
and State University, 1968. 

3. Youngdahl, C. K. : Influence of Pulse Shape on the Final Plastic Deforma- 

tion of a Circular Plate. Int. J. Solids Structures, Vol. 7, 1971, 
pp. 1127-1142. 

4. Florence, A. L.: Annular Plate Under a Transverse Line Impulse. AIAA 

Journal, Vol. 3, No. 9, 1965, pp. 1726-1732. 

5. Sneddon, I. E.: Fourier Transforms. McGraw-Hill Book Co., Inc., 1951, 

pp. 145-147. 

6. Madden, R. : Ballistic Limit of Double-Walled Meteoroid Bumper Systems. 

NASA TN D-3916, 1967. 

7. Thomson, R. G. : Plastic Behavior of Circular Plates Under Transverse 

Impulse Loadings of Gaussian Distribution. NASA TR R-279, 1968. 

8. Perzyna, P.: The Constitutive Equations for Work-Hardening and Rate 

Sensitive Plastic Materials. Proc. Vib. Problems, Vol. 4, 1963, 
pp. 281-290; Bull. Acad. Polon. Sci., Ser. Sci. Tech., Vol. 12, 1964, 
pp. 199-206. 

•• 

9. Hohenemser, K. ; and Prager, W. : Uber die Ansatze der Mechanik isotroper 

Konitinua. ZAMM, Vol. 12, 1932, pp. 216-226. 


611 



10. Prager, W. : Mechanique des solides isotropes an dela du domaine 

elastique. Memorial Sci. Math. 87, Paris, 1937, Eq. (47) on p. 27. 

11. Perzyna, P. : Fundamental Problems in Viscoplasticity. Advances in 

Applied Mechanics. Vol. 9, 1966, pp. 202-243. 

12. Wierzbicki, T.: Impulsive Loading of Rigid Viscoplastic Plates. 

Int. J. Solids Structures, Vol. 3, 1967, pp. 635-647. 

13. Florence, A. L. : Circular Plate Under a Uniformly Distributed Impulse. 

Int. J. Solids Structures, Vol. 2, 1966, pp. 37-47. 

14. Wierzbicki, T.: A Method of Approximate Solution of Boundary Value 

Problems for Rigid, Viscoplastic Structures. Acta Mechanica, Vol. 3, 
No. 1, 1967, pp. 56-66.. 

15. Hayduk, R. J. ; and Thomson, R. G. : Static Load-Carrying Capacity of 

Circular Rigid-Plastic Plates Under Gaussian Distribution of Pressure. 
Int. J. of Solids and Structures, Vol. 12, 1976, pp. 555-558. 

16. Wang, A. J.: The Permanent Deflection of a Plastic Plate Under Blast 

Loading. J. Appl. Mech. , Vol, 22, No. 3, 1955, pp. 375-376. 


612 




Figure 1.- Representation of pro- 
portional loading in deviatoric 
space. 



2 


Figure 2.- Simply supported circular 
plate with Gaussian distribution 
of initial velocity. 



Figure 3.- Element of the circular 
plate with internal forces and 
moments . 



0 .05 .10 .15 .20 .25 .30 


t, 10 3 sec 

Figure 4.- Representative time 
history of plate central 
velocity for the Gaussian 
ideal impulse loading. 
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(b) a - 1 x 10 sec. 

Figure 5.- Final plate profiles for 
various values of the ideal 
impulse shape parameter, 3* 




Figure 6.- Comparison of approximate 
expression (eq. (22)) for motion* 
to cease, tf , and points deter- 
mined from the complete equations 
for the ideal impulse loading. 



Figure 7.- Comparison of approximate 
expression (eq. (25)) for the 
final central deflection and points 
determined from the complete equa- 
tions for the ideal impulse loading. 



